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Abstract 

Background: Genomic imprinting is the epigenetic marking of genes that results in parent-of-origin monoallelic 
expression. Most imprinted domains are associated with differentially DNA methylated regions (DMRs) that originate 
in the gametes, and are maintained in somatic tissues after fertilization. This allelic methylation profile is associated 
with a plethora of histone tail modifications that orchestrates higher order chromatin interactions. The mouse 
chromosome 15 imprinted cluster contains multiple brain-specific maternally expressed transcripts including Ago2, 
Chracl, Trappc9 and Kcnk9 and a paternally expressed gene, Pegl3. The promoter of Pegl3 is methylated on the 
maternal allele and is the sole DMR within the locus. To determine the extent of imprinting within the human 
orthologous region on chromosome 8q24, a region associated with autosomal recessive intellectual disability, 
Birk-Barel mental retardation and dysmorphism syndrome, we have undertaken a systematic analysis of allelic 
expression and DNA methylation of genes mapping within an approximately 2 Mb region around TRAPPC9. 

Results: Utilizing allele-specific RT-PCR, bisulphite sequencing, chromatin immunoprecipitation and chromosome 
conformation capture (3C) we show the reciprocal expression of the novel, paternally expressed, PEG 13 non-coding 
RNA and maternally expressed KCNK9 genes in brain, and the biallelic expression of flanking transcripts in a range 
of tissues. We identify a tandem-repeat region overlapping the PEG! 3 transcript that is methylated on the maternal 
allele, which binds CTCF-cohesin in chromatin immunoprecipitation experiments and possesses enhancer-blocker 
activity. Using 3C, we identify mutually exclusive approximately 58 and 500 kb chromatin loops in adult frontal 
cortex between a novel brain-specific enhancer, marked by H3K4me1 and H3K27ac, with the KGNK9 and PEG 13 
promoters which we propose regulates brain-specific expression. 

Conclusions: We have characterised the molecular mechanism responsible for reciprocal allelic expression of the 
PEGU and KCNK9 transcripts. Therefore, our observations may have important implications for identifying the cause 
of intellectual disabilities associated with the 8q24 locus. 
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Background 

Genomic imprinting is the epigenetic marking of a sub- 
set of genes that results in parent-of-origin monoallelic 
expression. The regulation of imprinting is complex, 
involving interplay between many different epigenetic 
mechanisms, including DNA methylation, histone tail 
modifications and non-coding RNA (ncRNAs) [1]. Most 
imprinted genes are associated with a region of differen- 
tial DNA methylation that is acquired in the gametes, 
and maintained in somatic tissues after fertilization by 
the UHRF1-DNMT1 complex [2-4]. 

The majority of imprinted differentially DNA methyl- 
ated regions (DMRs) acquire their methylation from 
oocytes, as only a few examples of sperm-derived methy- 
lation at DMRs are known [5]. Maternally methylated 
DMRs generally act as promoters, associated in some 
cases with long ncRNAs. In some cases these ncRNAs 
confer silencing of neighbouring genes in cis through re- 
cruitment of histone remodelling complexes [6,7]. Some 
of the intergenic, paternally methylated, DMRs act as 
methylation-sensitive insulators recruiting CTCF [8], in- 
fluencing higher-order chromatin folding [8]. 

Aberrant imprinting in the brain is known to be asso- 
ciated with severe developmental disorders such as 
Angelman and Prader-Willi syndromes, and other be- 
havioural phenotypes including autism-spectrum dis- 
order (ASD), mental retardation (MR) and psychosis 
[9,10]. Recently, several genome-wide scans have identi- 
fied susceptibility alleles for MR and ASD mapping to 
8q22-24 [10]. This approach was used to identify mater- 
nally inherited missense mutations of the imprinted two 
pore-domain potassium channel (K2P) gene, KCNK9 
(also known as TASKS) in Birk-Barel mental retardation 
syndrome [11]. This gene is located adjacent to 
TRAPPC9 (also known as NIBP), a gene disrupted by 
both coding mutations and deletions in non-syndromic, 
autosomal-recessive mental retardation [12-15]. The 
orthologous region on mouse chromosome 15 shows 
parent-of-origin allelic expression in mouse brain, and 
includes the maternally expressed Kcnk9 [16], Trappc9, 
Chracl and Ago2 (also known as Eif2c2) genes, clustered 
around the paternal Pegl3 ncRNA (Figure 1A). 

In the present study, we characterise the human 8q24 
locus, and show that only PEG 13 and KCNK9 are 
imprinted (Figure IB). These transcripts are reciprocally 
expressed, with paternal expression of PEG13 and ma- 
ternal expression of KCNK9. The PEG 13 transcript is 
embedded within a maternal methylated region, with 
the unmethylated allele enriched for CTCF-cohesin, 
which acts as an enhancer-blocker in vivo. Lastly, we 
show chromatin loops between an enhancer region, 
with a canonical enhancer histone modification signa- 
ture in brain tissues, and the KCNK9 and PEG13 
promoters. 



Results 

Imprinting within human chromosome 8q24 

A 700 kb region on mouse chromosome 15 was recently 
identified as a novel imprinted domain using ultra- 
sensitive RNA-seq technology [17]. The cluster of four 
maternally expressed genes, as well as numerous 
expressed sequence tags (ESTs) is located around the pa- 
ternally expressed 4.7 kb ncRNA, Pegl3, within intron 
17 of the Trappc9 gene. The Trappc9 gene itself is 
subject to isoform-specific imprinting, with a single 
truncated isoform being paternally expressed, however, 
all the full-length transcripts are maternally derived 
(Figure 1A) [17]. The promoter of Pegl3 is a CpG island 
(CpG44 in Figure 1A) that is methylated on the maternal 
allele, being established in developing oocytes [18,19]. 

To assess if imprinting is conserved within the human 
orthologous region, we assessed the allelic expression for 
COL22A1, KCNK9, TRAPPC9, CHRAC1 and AG02 
(Figure IB). We identified transcribed single nucleotide 
polymorphisms (SNPs) that would allow for allelic dis- 
crimination, and imprinting analysis was performed in 
adult leucocytes, brain and term placentae. We confirm 
that KCNK9 is monoallelically expressed in the majority 
of brain regions, from the maternal allele (n = 2; includ- 
ing whole brain, hippocampus, cerebellum, vermis, par- 
ietal lobe and entorhinal cortex). RT-PCR followed by 
pyrosequencing confirmed the Sanger sequencing results 
and revealed that only residual expression arises from 
the repressed allele (Additional file 1: Figure SI). All 
other transcripts were expressed biallelically in multiple 
tissues (Figure ID; Additional file 2: Table S2). 

The human PEG13 ncRNA is comprised of numerous 
expressed tandem-repeats 

We identified numerous ESTs (reference AX748239), the 
majority derived from the cerebellar and amygdalar re- 
gions of the brain, overlapping a large CpG island within 
intron 17 of the human TRAPPC9 gene (CpG210 in 
Figure IB), which is in the orthologous intron as the 
mouse Pegl3 gene. Sequence analysis revealed no simi- 
larity with mouse Pegl3, despite their similar location. 
The region is made up of numerous degenerative simple 
tandem repeats, including 4.7 copies of a 96 bp repeat 
and 43.5 copies of a 15 bp repeat. Despite numerous at- 
tempts, we failed to detect a signal by northern blot ana- 
lysis, however this region is associated approximately 6 
kb of contiguous overlapping RNA-seq reads in brain 
derived RNA (Burge lab brain and ENCODE RNA-seq 
track UCSC genome browser GRCh37/hgl9) suggesting 
that expression across this repetitive region produces a 
single transcript (Figure 1C). Using qRT-PCR we were 
able to determine that the expression was most abun- 
dant in brain (Additional file 1: Figure SI) and, in sam- 
ples heterozygous for rs4289794, was monoallelically 
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Figure 1 Schematic representation of the Peg13/PEG13 imprinted domain on distal mouse chromosome 15 and human chromosome 
8q24. (A) Map of the Troppc9/Pegl3 locus located on mouse chromosome 15, showing the location of the various imprinted transcripts 
(red transcripts are maternally expressed, blue are paternally expressed and black are expressed from both parental alleles; arrows represent 
direction of transcription) and CpG islands (green bars). (B) Schematic of the human TRAPPC9/PEG1 3 region on chromosome 8. (C) RNA-seq 
profiles throughout the human domain in brain and lymphoid tissue. The RNA peaks corresponding to PEG! 3 and KCNK9 are indicated in the 
brain track. The RNA density plots represent relative read density (normalised counts per million mapped reads as defined by the University 
of California Santa Cruz (USCS) sequence browser). (D) The allelic expression for the genes in the human cluster. The sequence traces for 
heterozygous DNA samples from brain, leucocytes and term placentae are shown, as well as the resulting RT-PCR. 



expressed (n = 10; including cerebellum, vermis, hippo- 
campus, parietal lobe and occipital cortex). Pyrosequenc- 
ing assay confirmed these observations (Additional file 
1: Figure SI). One sample was accompanied by inform- 
ative parental DNA samples and so the paternal origin 
of the expressed allele was ascertained (Figure ID). 



The CpG island overlapping PEG13 is a maternally 
methylated DMR 

It has previously been proposed that imprinting of this 
cluster in mouse is controlled by the maternally methyl- 
ated germline DMR encompassing the promoter CpG 
island associated with the Pegl3 [18]. We analysed the 
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DNA methylation status of all human promoter associ- 
ated CpG islands throughout the domain in human leu- 
cocytes and brain tissues (grey and white matter isolated 
from the cerebral cortex, specifically the dorsolateral pre- 
frontal cortex; Brodmann area 9) using a combination whole 
genome bisulphite sequencing (WGBS) and ChlP-seq for 
meDIP and H3K4me3, a histone modification refractory to 
DNA methylation. This high-throughput analysis revealed 
that the KCNK9, AG02, CHRAC1 and TRAPPC9 promoter 
CpG islands are unmethylated in both tissue types and are 
associated with abundant H3K4me3 enrichment (Figure 2A). 
In addition, we observe that the CpG island overlapping 
PEG13 is approximately 50% methylated in brain and 
lymphocyte WBGS datasets and is associated with co- 
enrichment of H3K4me3 and meDIP, consistent with differ- 
ential active and repressive chromatin states on homologous 
chromosomes (Figure 2A). We subsequently analysed the 
methylation status of this region using bisulphite PCR and 
direct sequencing. This revealed that the PEG13 region 
flanking the rs4455807 SNP exhibits allelic methylation in 
brain (whole fetal and adult cerebellum, hippocampus, 
frontal and occipital cortex), leucocytes and placenta- 
derived DNA (Figure 2B). In informative cases, the DMR 
was methylated on the maternal allele. In seven brain sam- 
ples heterozygous for this polymorphism, we observed ro- 
bust allelic methylation and monoallelic expression of the 
unmethylated allele, consistent with paternal expression 
(Figure 2C). One of these samples was directly shown to be 
paternally expressed when parental DNA samples were in- 
terrogated. Coherent with this observation, this interval is 
devoid of methylation in sperm derived DNA (Figure 2B). 

The paternal allele of the PEG 73-DMR binds CTCF-cohesin 

Using a panel of DNA samples derived from various 
brain regions (parietal lobe, occipital lobe, cerebellum, 
frontal cortex, temporal gyrus, hippocampus and vermis) 
we confirm the unmethylated status of the promoter of 
KCNK9, which agrees with the observations obtained 
from WGBS (data not shown). This suggests that the 
imprinting of this gene is dependent on other ds-acting 
regulatory elements. The maternal expression of KCNK9 
and the maternal methylation of the PEG13-DMR are 
consistent with the action of an enhancer-blocker mech- 
anism utilizing CTCF, similar to that which regulates the 
reciprocal imprinting of Igf2/H19 [20]. Several strong ca- 
nonical two-part CTCF (motif 1 + 2) binding sites within 
the PEG13-DMR were revealed using an in silico analysis 
using the published ChlP-sequence data (data not 
shown) [21,22]. To confirm in vitro binding, we investi- 
gated CTCF binding in published ChlP-seq datasets. We 
observed CTCF enrichment at the PEG13-DMR in brain 
and lymphocyte experiments and confirmed these obser- 
vations utilizing ChIP on normal lymphoblast cells 
(region 2 of Figure 3). The efficiency of the ChIP was 



confirmed by parallel analysis of the 7/79-ICR, with ma- 
ternal enrichment observed for a PCR amplicon incorp- 
orating CTCF sites 4 to 7 (Figure 3C) [23]. Subsequent 
analysis revealed precipitation of the PEG13-DMR. Un- 
fortunately, the cell lines were not informative for any 
polymorphisms that would allow discrimination of par- 
ental alleles. To circumvent this, we performed bisulphite 
pyrosequencing for the PEG13-DMR on CTCF immuno- 
precipitated DNA. This revealed that, similar to the H19- 
ICR, CTCF preferentially binds to the unmethylated allele 
(Figure 2D). Recently, cohesin has been shown to play a 
critical role in maintaining CTCF higher-order chromatin 
confirmation at the H19-IGF2 loci [8,24]. To determine if 
CTCF and cohesin co-localize at the PEG13-DMR, we 
performed ChIP using antisera against the cohesin sub- 
units RAD21 and SMC3. In lymphoblast cell lines, the 
cohesin subunits clearly precipitated at the PEG13-DMR 
(Figure 3B), suggesting that a CTCF-cohesin complex may 
orchestrate higher order chromatin loops. 

The PEG 73-DMR has enhancer-blocking 

To test the PEG13-DMR for enhancer-blocking activity 
we used a previously described enhancer-blocking assay 
[25,26]. We studied the ability of a 480 bp PEG13-DMR 
fragment containing the canonical CTCF motifs to inter- 
fere with the activity of a heterologous enhancer/pro- 
moter interaction in HEK293 cells. The PEG13-DMR 
construct showed consistent enhancer-blocking activity, 
represented as fold reduction of activity compared to a 
vector lacking an insert (Figure 4). The PEG13-DMR 
fragment showed a stronger enhancer-blocking activity 
than the II/III control fragment containing the (3-globin 
HS4 core' enhancer-blocker, but not as much as the 
5HS4 insulator element. 

CTCF mediated allele-specific chromatin looping occurs 
between a shared enhancer and the KCNK9 and PEG13 
promoters 

Recent studies have revealed that CTCF occupancy corre- 
lates with domain boundaries, as well as mediating intra- 
and interchromosomal contacts [27]. In an attempt to 
understand the mechanism regulating reciprocal brain- 
specific expression of KCNK9 and PEG13, we interrogated 
publically available genome-wide ChlA-PET from MCF-7 
cells and CTCF ChlP-seq datasets. We identified add- 
itional strong CTCF enrichment at the KCNK9 promoter 
and a region within intron 17 of TRAPPC9 (regions 1 and 
3 of Figures 3 and 5C, respectively), which has an enhan- 
cer chromatin signature with the co-enrichment of 
H3K4mel, H3K27ac and p300 in various brain ChlP-seq 
datasets, including the frontal cortex [28,29] . Interestingly, 
in lymphocytes, which do not express KCNK9 or PEG13, 
this region is not associated with either histone modifica- 
tion, suggesting that that it is a brain-specific enhancer 
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(See figure on previous page.) 

Figure 2 DNA methylation profiling of gene promoter of genes flanking TRAPPC9/PEG13. (A) CpG methylation profiling in lymphocytes 
and brain samples as defined by WGBS. The analyses were restricted to the intervals overlapping the promoters/CpG islands and are associated 
with the expected H3K4me3 and meDIP signatures in CD4 lymphocytes cells and fetal brain. The two shades of grey peaks in the meDIP 
and ChlP-seq panels represent two independent biological replicates compared to input (black peaks) with the y-axis showing the number 
of ChlP-seq reads. The precise location of the PEG13-DMR, as defined by approximately 50% methylation and co-enrichment of H3K4me3 and 
meDIP is indicated by the red bracket. (B) The methylation status of the PEG13-DhAR was confirmed using standard bisulphite PCR on DNA 
samples heterozygous for SNP rs4455807. (C) Sequence traces for two brain samples revealing the genotypes of the PEG13-DMR methylated 
allele (resistant to Hha\ digestion) and expression from the unmethylated allele. 



(Figure 5C). Subsequent ChIP analysis confirmed that 
both of these regions were strongly enriched for both 
CTCF and the cohesin subunits RAD21 and SMC3 
(Figure 3B), both of which had insulator function in our 
enhancer-blocking assay (Additional file 3: Figure S2). The 
ChlA-PET dataset revealed that these CTCF-cohesin re- 
gions, separated by approximately 500 kb, physically inter- 
act (Figure 5A). 

To confirm this physical interaction, we performed 
chromatin conformation capture experiments (3C) to 
identify potential chromatin folding. 3C-qPCR assays 
were performed on cerebellar samples and interaction 
frequencies were determined between a constant Hindlll 
site located within the unmethylated CTCF-cohesin 
binding site within the KCNK9 promoter and other Hin- 
dlll sites throughout the locus. We identified strong in- 
teractions between the KCNK9 promoter constant 
fragment with the PEG13-DMR and the CTCF-cohesin 
site in the enhancer region located within intron 17 of 
TRAPPC9 (Figure 5B). Direct sequencing of the chimer- 
ical products from brain revealed that the appropriate 
chimeric products result from the 3C ligations (data not 
shown). We subsequently used a PCR primer within the 
enhancer region as the constant, to not only confirm the 
physical interaction with the KCNK9 promoter, but to 
also determine whether this region juxtaposes the 
PEG 13 promoter. Using this second constant fragment, 
we observed that the enhancer region physically interacts 
with both the KCNK9 and PEG13 in human cerebellar 
samples (Figure 5B). Unfortunately, no informative SNPs 
were identified within the vicinity of the Hindlll sites asso- 
ciated with the enhancer that would allow us to determine 
the parental origin of the resulting chemical products. 
These results suggest that the brain-specific reciprocal 
imprinted expression of PEG 13 and KCNK9 is dependent 
upon the active enhancer configuration and higher-order 
chromatin looping (Figure 5D). 

Discussion 

We have characterised a new imprinted cluster on hu- 
man chromosome 8q24, which harbours two brain- 
specific genes. We show that the promoter of the novel, 
paternally expressed, non-coding PEG13 transcript arises 
from a maternally methylated DMR. This DMR binds a 



CTCF-cohesin complex that acts as a methylation- 
sensitive enhancer-blocker on the unmethylated paternal 
allele. Furthermore, we have utilized publically available 
Chi A- PET datasets to identify an enhancer region that 
was previously unknown to be associated with gene ex- 
pression within the locus. The ChlA-PET technique 
combines both chromatin immunoprecipitation with 
3C-type analyses for the direct analysis of chromatin in- 
teractions exclusively formed between sites bound by a 
given chromatin interacting protein, in our case CTCF, 
however it does not address whether the identified loops 
depend on the protein of interest. We confirm the 
Chi A- PET interactions using standard 3C assays in 
brain. Our data suggest that the brain-specific expression 
of the imprinted transcripts is directly due to the acquisi- 
tion of an active enhancer chromatin signature. The gen- 
omic regions immediately flanking the 3C/ChIA-PET 
interacting region within intron 17 of TRAPPC9 is marked 
in neuronal tissues by H3K4mel, a histone mark associ- 
ated with many cell-type specific enhancers, abundant 
H3K27ac, with the acetyl group presumably deposited by 
the P300 acetyltransferase [30]. 

DNA methylation composition in brain is highly dy- 
namic, being present in several forms. An oxidative 
derivative of 5-methylcytosine (5mC) called 5-hydroxy- 
methylcytosine (5-hmC) accumulates in adult brain 
[31,32]. In addition, methylation in non-CpG context 
(mCH, where H is A, C or T) is also more abundant in 
frontal cortex of both adult mouse and humans [32]. 
Previous reports have shown that the adult brain con- 
tains high levels of 5-hmC, accounting for approximately 
40% of CG methylation in cerebellar Purkinje cells [31]. 
This modification is especially enriched in multi-exonic, 
highly expressed brain genes [33]. Methodologies based 
on standard bisulphite-sequencing techniques cannot 
distinguish between abundant 5mC and 5-hmC and 
sites, so the methylation profiles we identified by WGBS 
or bisulphite PCR will represent the sum of both modifi- 
cations. However, our results using meDIP using an anti- 
body directed against 5mC strongly suggest that the 
methylation we observe is not 5-hmC. This observation 
at the brain-specific PEG13 transcript is consistent with 
intronless, single exon genes having significantly less, or 
no, 5-hmC when compared to spliced transcripts [34]. 
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Figure 3 ChIP analysis for CTCF and the cohesion subunits RAD21 and SMC3. (A) ChlP-Seq data analysis in lymphocytes and cerebellum 
reveals the location of three ubiquitous CTCF binding sites. The positions of the ChIP PCRs are indicated. (B) qPCR performed on CTCF, RAD21 
and SMC3 ChIP material in normal lymphoblastoid cells at the intervals identified by ChlP-Seq. Graphs are represented as % of precipitation 
relative to input chromatin (mean values ± SEM). (C) Sequence traces showing monoallelic precipitation of CTCF and cohesion subunits at the 
control H79-ICR. (D) The methylation levels of at the H19-ICR and PEG13-DMR in CTCF input and ChIP material as determined by bisulphite PCR 
followed by pyrosequencing. 



In addition, a previous genome-wide study has identified 
'diffuse' parent-of-origin dependent non-CG methylation 
sites overlapping the gene bodies of the mouse Ago2 and 
Pegl3 genes [19]. Our analyses did not reveal any allelic 
methylation outside the CpG context at the PEG 13 or 
AG02 gene in WGBS datasets from human grey and 



white brain matter or by direct bisulphite PCR in adult 
frontal cortex (Additional file 4: Figure S3) [19]. This 
evolutionary discrepancy may partially explain the lack 
of imprinting of the AG02 gene in humans. 

Changes in the accumulation of DNA methylation in 
neurons throughout early childhood has been implicated 
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in learning and memory, as well as cognitive function 
[32,35]. Some Angelman syndrome cases, as well as aut- 
ism and ASD, are associated with coding mutations and 
methylation defects at the imprinted SNRPN locus 
[9,10]. Interestingly, similar to the long-range looping we 
that we have identified for KCNK9, the PWS/AS-IC in- 
teracts with the CHRNA7 gene, a distance of more than 
5 Mb [36,37]. This higher-order chromatin interaction is 
associated with MeCP2, suggesting that this protein, like 
CTCF, is capable of orchestrating long-range chromatin 
organisation in the brain [37]. 



Conclusions 

We identified a maternally methylated region that binds 
CTCF-cohesin, possesses enhancer- blocker activity, 
which we hypothesize dictates mutually exclusive chroma- 
tin looping between a novel enhancer region and the pro- 
moters of the reciprocally imprinted PEG 13 and KCNK9 
transcripts. Our findings are of particular relevance, not 
only to those researchers working on the aetiology of in- 
tellectual disabilities, but also to those interested in tissue- 
specific transcriptional regulation, since we show that the 
essential chromatin scaffold of the domain is present in all 
tissues but that a brain-specific enhancer is the only regu- 
latory element dictating tissue specificity. It will be inter- 
esting to determine whether the PEG 13 ncRNA functional 
contributes to the imprinting of the KCNK9 gene and 
whether these two imprinted transcripts are subject to 
epigenetic deregulation of the ds-acting sequences that 
could silence or cause loss-of-imprinting in non- 
syndromic forms of intellectual disability. 



Methods 

Human tissues 

A cohort of 48 placental samples with corresponding 
cord blood was collected at Hospital Sant Joan de Deu 
(Barcelona, Spain). DNA and RNA extraction and cDNA 
synthesis were carried out as previously described [38]. 
Normal peripheral blood was collected from adult vol- 
unteers aged between 19 and 60 years of age. A total of 
40 normal adult brain samples were obtained from 
BrainNet Europe/Barcelona Brain Bank. The dissection 
of individual brain regions (hippocampus, cerebellum, 
vermis, parietal lobe, entorhinal and occipital cortices) 
was performed by an experienced pathologist on ca- 
davers within 14 hours of death. The human RNA panel 
was purchased from Clontech (Human Total RNA mas- 
ter Panel II). Ethical approval for collecting cord blood, 
placental biopsies, brain samples and adult blood sam- 
ples was granted by the ethical committee of Hospital 
Sant Joan de Deu Ethics Committee (study number 35/ 
07) and IDIBELL (PR006/08 and PR048/13). 



Cell lines 

Control lymphoblastoid cell lines were established by 
EBV transformation of peripheral blood cells and propa- 
gated as previously described [7]. Prior to ChIP, the 
lymphoblastoid methylation signature throughout the 
8q24-imprinted domain was compared to leucocytes to 
ensure that the transformation process had not altered 
the epigenetic profile. The HEK293 cell line was used 
for the enhancer-blocking assay and was grown in 
DMEM supplemented with 10% FCS and antibiotics. 
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Figure 5 Chromatin interactions with CTCF sites upstream of the KCNK9 promoter, the PEG 73-DMR and a novel brain-specific enhancer. 

(A) A schematic map of the TRAPPC9 and KCNK9 genes and annotated ChlA-PET data showing CTCF interactions. The red star represents the 
location of the brain-specific enhancer. (B) The 3C looping profile in brain using a constant primer 5' to the KCNK9 promoter (upper graph) or 
within the enhancer located within intron 17 of TRAPPC9 (lower graph). (C) Sequence track for H3K4me1 and H3K27ac for a 30 kb region 
containing the TRAPPC9 intron 17 ChlA-PET site. Tracks are shown in expressing in brain and non-neuronal lymphocyte samples. The two shades 
of grey peaks in the meDIP and ChlP-seq panels represent two independent biological replicates compared to input (black peaks) with the y-axis 
showing the number of ChlP-seq reads. The red star depicts the location of the proposed enhancer region. (D) A proposed model of higher- 
order chromatin looping within the 8q24 domain, showing the special organisation of the allele-specific PEG13 enhancer-blocker in relation to 
the gene promoters, CTCF binding sites and H3K4me1/H3K27ac enhancer region. (•) methylated cytosine, (O) unmethylated cytosines and green 
block representing the H3K4me1/H3K27ac region associated with the brain-specific enhancer (e). 
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Epigenetic bioinformatics analysis 

RNA-seq datasets were downloaded from the NCBI 
GEO repository (GSM325476 and GSM 325483). Gen- 
omic data for H3K4me3 and corresponding meDIP in 
brain and non-neuronal samples used were retrieved 
from GSM54305, GSM613913, GSM772916, GSM772 
948, GSM772836, GSM669615, GSM669614, GSM806 
948, GSM806935 and GSM806943. The brain and leuco- 
cyte whole genome bisulphite sequencing datasets were 
from Lister et al [32] (GSE47966, GSE46698 and 
GSE31263). The human CTCF, H3K4mel, H3K27ac and 
P300 ChlP-seq and ChlA-PET datasets were those with 
visual tracks on the UCSC genome browser (GRCh37/ 
hgl9) and raw reads were retrieved from GSM1022661, 
GSM1022662, GSM1022673, GSM651541, GSM772884, 
GSM772924, GSM772930, GSM772881, GSM772963, 
GSM772997, GSM670036, GSM772992, GSM772991, 
GSM670010, GSM77295, GSM773007, GSM670033, 
GSM772990, GSM670003, GSM772985, GSM772985 
and GSM970215. Genomic mapping of all reads and 
data analysis was done with in-house scripts using the R 
statistical package. 

Allelic expression analysis 

The placental and brain tissue genotypes of the expressed 
genes were obtained by PCR and direct sequencing. Se- 
quences were interrogated using Sequencher v4.6 (Gene 
Codes Corporation, MI, Ann Arbor, USA) to distinguish 
heterozygous and homozygous samples. Heterozygous 
sample sets were analysed for allelic expression using RT- 
PCR that incorporated the polymorphisms in the final 
PCR product (Additional file 5: Table SI). The resulting 
RT-PCR amplicons were sequenced in both directions. 
The amplification cycle numbers for each transcript were 
determined to be within the exponential phase of the 
PCR, which varied for each gene, but was between 32 and 
38 cycles. 

DNA-methylation analysis 

Approximately 2 ug DNA was subjected to sodium 
bisulphite treatment and purified using the EZ GOLD 
methylation kit (ZYMO, Orange, CA, USA) and was 
used for bisulphite PCR analysis. Bisulphite PCR primers 
for each region were used with Hotstar Taq polymerase 
(Qiagen, Crawley, UK) at 40 cycles and the resulting PCR 
product cloned into pGEM-T easy vector (Promega, 
Madrid, Spain) for subsequent sequencing (Additional file 
5: Table SI). 

Real-time RT-PCR 

All PCR amplifications were run in triplicate on either 
Applied Biosystems 7500 or 7900 Fast real-time PCR 
machines (Applied Biosystems, Life Technologies, Grand 
Island, NY, USA) following the manufacturers' protocol. 



All primers (Additional file 5: Table SI) were optimized 
using SYBR Green (Applied Biosystems, Life Technolo- 
gies, Grand Island, NY, USA) and melt curve analysis to 
ensure that amplicons were specific and free of primer- 
dimer products. Thermal cycle parameters included Taq 
polymerase activation at 95°C for ten minutes for one 
cycle, repetitive denaturation at 95°C for 15 seconds, 
and annealing at 60°C for one minute for 40 cycles. All 
resulting triplicate cycle threshold (Ct) values had to be 
within one Ct of each other for the sample to be in- 
cluded. The quantitative values for each triplicate were 
determined as a ratio of LI 9 expression, which was mea- 
sured in the same sample, and the triplicate mean used 
to provide relative expression values. 

Pyrosequencing 

Pyrosequencing was used as an accurate method of 
quantifying allelic expression of PEG 13 and KCNK9 in 
heterozygous brain samples. Standard RT-PCR was used 
to amplify across SNPs with the exception that the re- 
verse primers were biotinylated. The entire biotinylated 
RT-PCR product (diluted to 40 ul) was mixed with 38 ul 
of binding buffer and 2 ul (10 mg/ml) streptavidin- 
coated polystyrene beads. Bead-amplicon complexes 
were captured on a vacuum prep tool (Qiagen, Crawley, 
UK) and the PCR products denatured using 0.2 M 
NaOH. The denatured DNA was resuspended in 40 
pmol of sequencing primer dissolved in 12 ul water and 
primers annealing was achieved by heating the sample 
to 80°C for two minutes before cooling to room 
temperature. For sequencing, forward primers were de- 
signed to the complementary strand (Additional file 5: 
Table SI). The pyrosequencing reaction was carried out 
on a PyroMark Q96 instrument (Qiagen, Crawley, UK). 
The peak heights were determined using the pyrose- 
quencing commercial software. For methylation pyrose- 
quencing of the PEG13-DMR, the same protocol was 
followed with the exception that the initial template was 
amplified from a bisulphite converted template and the in- 
terrogated sites were C/T variants at CpG dinucleotides. 

Chromatin immunoprecipitation (ChIP) 

Chromatin from about 80 million cells was aliquoted 
into 100 mg batches and used for each immunoprecipi- 
tation reaction with Protein A Agarose/Salmon Sperm 
DNA (Millipore, 16-157, Billerica, MA, USA) and spe- 
cific antibody. The antibody against CTCF (07-729) was 
obtained from Millipore (Billerica, MA, USA) and the 
antibodies against RAD21 (AB992) and SMC3 (AB9263) 
were obtained from Abeam (Cambridge, MA, USA). For 
each ChIP, a fraction of the input chromatin (1%) was 
also processed for DNA purification and a mock immu- 
noprecipitation with a neutral, unrelated IgG antiserum 
was carried out in parallel [24]. 
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Levels of immunoprecipitated chromatin at specific re- 
gions were determined by qPCR with an Applied Biosys- 
terns 7900 Fast real-time PCR machine, using SYBR 
Green (Applied Biosystems, Life Technologies, Grand Is- 
land, NY, USA). Each PCR was run in triplicate and pro- 
tein binding was quantified as a percentage of total 
input material 

Enhancer-blocker/lnsulator assay 

We used an enhancer-blocking assay to address the insula- 
tor activity of a DNA fragments containing the PEG13- 
DMR, KCNK9 promoter and the ChlA-PET enhancer 
region utilizing the pELuc plasmid system [26]. The DNA 
fragments to be tested were cloned between the cyto- 
megalovirus (CMV) enhancer and the promoter (Xhol) or 
upstream of the CMV enhancer (Pstl). The assay was per- 
formed by transfecting the constructs into HEK293 cells. 
The resulting data are presented as fold-enhancer-blocking 
activity normalized to the values achieved by the basal 
pELuc vector. As positive controls, we used the 5HS4 (1.2 
kb chicken p-globin insulator) or the II/III (5'HS4 core' 
p-globin) insulator element [25]. For a negative control, we 
used II/III fragment with mutated CTCF binding sequence. 

3C analysis 

The chromatin conformation capture (3C) protocol was 
performed as previously described [39] with minor 
amendments. Approximately 100 mg of snap frozen 
brain tissue was reduced to powder with a pestle and 
mortar under liquid nitrogen. The pulverized brain sam- 
ples were cross -linked with 1% formaldehyde for ten mi- 
nutes and the reaction was blocked by adding glycine to 
a final concentration of 0.125 M. Subsequently, Hindlll 
was used to digest 1 x 10 7 formaldehyde cross -linked nu- 
clei (overnight digestion, 1,200U, New England Biolabs, 
Ipswich, MA, USA ). The efficiency of the restricted en- 
zyme digestion was assessed by qPCR across each re- 
striction site, comparing digested and undigested 
chromatin fractions. Only chromatin with digestion effi- 
ciency above 80% was used. Subsequently, the DNA was 
ligated overnight in a 500 \A reaction volume using 
1,950 units of T4 ligase (Fermentas- Fisher Scientific, 
Madrid, Spain). DNA was decross -linked by incubating 
overnight at 65°C and purified using phenol/chloroform 
extraction. This DNA was used for RT-PCR (LightCy- 
cler, Roche Applied Science, Barcelona, Spain) to deter- 
mine the frequency of interactions, using constant 
primers either in the KCNK9 promoter or the enhancer 
(Additional file 5: Table SI). Primer efficiency and basal 
interaction frequencies were determined using digested 
and ligated bacterial artificial chromosome (BAC) DNAs 
(human RPII 1069118 and RPII 431 L20) as described by 
Braem and co-workers [39]. 



Additional files 

Additional file 1: Figure SI. (A) The expression of TRAPPC9, PEGU and 
KCNK9 in a panel of human tissues as determined by qRT-PCR. All values 
are relative to the housekeeping gene RPL19. (B) The confirmation of 
KCNK9 and PEGU allelic expression using pyrosequencing. 

Additional file 2: Table S2. The number of heterozygous tissue 
samples used to determine allelic expression of novel imprinted 
transcripts. 

Additional file 3: Figure S2. Determining the insulator activity of 
additional CTCF ChlA-PET regions. The bars indicate the firefly luciferase 
expression relative to Renilla luciferase activity for constructs containing a 
370 bp fragment encompassing the CTCF within the KCNK9 promoter or 
a 510 bp fragment containing the CTCF adjacent to the brain-specific 
enhancer. As a control, the enhancer-blocking assay was also performed 
with inserts for the 5HS4 (1.2 kb insulator), the II/III ('core' 5'HS4) of 
(3-globin enhancer-blocker, as well as the II/III construct with mutated 
CTCF sites (II/III mut). The constructs are illustrated at the top of the 
figure. Data are presented as fold-enhancer-blocking activity normalised 
to the reference pELuc vector. The experiment represents the means of 
triplicate reading (±SD) with independent replicate experiments giving 
comparable results (data not shown). 

Additional file 4: Figure S3. The analysis of non-CG methylation within 
the gene bodies of PEGU and AG02. (A) The CpG methylation (black) 
and CH methylation (grey) were determined from WGBS for grey and 
white matter. Only background levels of non-CG methylation were 
observed, equating to 1.3% and 0.7% at PEGU and 1.1% and 0.6% at 
AG02 in grey and white matter respectively. (B) Bisulphite sequence in 
frontal cortex derived DNA samples. The large circles represent single 
CpG dinucleotides whereas smaller circles depict individual CH on the 
strand, (•) a methylated cytosine, (O) unmethylated cytosines. 

Additional file 5: Table SI. Primers used in this study. 
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